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Abstract 

We propose a multi-graviton theory with non-nearest-neighbor couplings in the the- 
ory space. The resulting four-dimensional discrete mass spectrum reflects the structure 
of a latticized extra dimension. For a plausible mass spectrum motivated by the dis- 
cretized Randall-Sundrum brane-world, the induced cosmological constant turns out 
to be positive and may serve as a quite simple model for the dark energy of our accel- 
erating universe. 
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1. Introduction. There is renewed interest in the study of multi-graviton theories [1]. 
Generally speaking, the corresponding action can be viewed as a sum of Pauli-Fierz 
terms while a consistent interaction of gravitons does not occur there. In a way, multi- 
graviton theories resemble higher-dimensional gravities with a discretized dimension. 
And such sort of discretized Kaluza-Klein theories are now under intense attention 
due to their primary importance in the realization of the dimensional deconstruction 
program [2, 3]. Moreover, multi-gravitons are also related with discretized brane- worlds 
[4]- 

In spite of the the absence of a consistent interaction among the gravitons, one can 
think on possible couplings in the theory space. In particular, in a recent paper [5] a 
multi-graviton theory with nearest-neighbour couplings in the theory space has been 
proposed. As a result, a discrete mass spectrum appears. The theory seems to be 
equivalent to Kaluza-Klein gravity with a discretized dimension. 

In this letter, we propose a simple model of a multi-graviton theory with non- 
nearest-neighbour couplings in the theory space. Depending on the choice for such 
couplings, a the different graviton mass spectrum is induced. For instance, our model 
may be seen as motivated by a discretized Randall-Sundrum (RS) brane-world [6] or 
by complicated boundary condition on a latticized dimension. It opens number of 
possibilities for the choice of mass spectrum. For example motivated by discretized RS 
brane-world model the small positive cosmological constant is induced as shown by the 
explicit one-loop calculation. 

2. A multi-graviton model. We consider N variables, (f> n , which may be identified 
with the fields on a lattice with N sites. The difference operator A is defined by 

JV-l 

A0„ = ^2 a k 4> n +k • (1) 

Here we assume 4> n +N = 4>m which may be regarded as a periodic boundary condition. 
In (1), the afc's are N constants. In order that A becomes a usual differentiation in a 
proper continuum limit, the following condition is imposed: 

N-l 



«fe = • ( 2 ) 



k=0 

The eigenvectors for A are given by 

0f= e*T, M = 0,l,2,3 (3) 
v jV 

and their corresponding eigenvalues, by 

N-l 

A^ = im M ^, ™ M = ^a n e^. (4) 



n=0 



Here m° = 0. We should also note that (f>^ satisfies the following properties, which 
may be identified with the conditions of normalization and completeness, respectively: 



N-l N-l 



<Pn 4>n =0 , l^<?n <Pn> = b nn< ■ (5) 



n=0 M=0 
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One can solve a n with respect to m M by 

. N-l . N-l 

1 E = i E - M e"^ • (6) 



M=0 M=0 

Then by choosing a n properly, we may obtain arbitrary spectrum of m M with m° = 0. 
We should note if a n is real 

m N - M = -(m M )*. (7) 

In what follows, we will assume both that a n is real and that the condition (7) is 
satisfied. For instance, if N = 2, we find that m 1 is pure imaginary and that 

ao = -9m 1 , a\ = 9m 1 . (8) 

Here 9m 1 is the imaginary part of m 1 , and we denote the real part of m by 3ftm. For 
N = 3, we have m 2 = — (m 1 )*, and 

ao = — -9m 1 , ai = - ^Sm 1 + \/3lftmi^ , 

a 2 = - (3m 1 - \/3Kmi) . (9) 

For N = 4, we have m 3 = — (m 1 )* and find m 2 is pure imaginary. Then we obtain 

a = \ (29m 1 - 9m 2 ) , ai = \ ( 2Kml + 9m2 ) > 

o 2 = \ (-29m 1 - 9m 2 ) , a 3 = - (-2Km 1 + 9m 2 ) . (10) 

One may apply the above arguments to the multi-graviton theory , extending the 
formalism in [5]. 

The Lagrangian of the massless spin-two field (graviton) h^ v is given by 

-Co = -\d x h^d x h^ + d x h\d u h^ - d^ u d u h + ^d x hd x h , (11) 

where h = W 1 ^. The Lagrangian of the massive graviton with mass m is given by 
introducing Stuckerberg fields and ip [7] 



2 

C m = C - ^- ( h^hT - h 2 ) - 2 (m# + &>tp) (d u V " 9h) 

~ (dpAv - d„Ap) (d^A v - d v A») . (12) 

We consider N copies of the graviton fields h niiv and Stuckerberg fields A nil and ip n . 
By replacing 

h^ v -> Kin, , \ -> —iA nil , -> y?„ , m -> -iA, (13) 

and adding all terms in the action (sum on n), we hereby propose a theory given by 
the following Lagrangian, which is a generalization of the one in [5]: 

N-l r 

-^QxK^d x K: + d x h x fijc - d^K^dvK + -d x h n d x h n 



£ = E 

n=0 

+ l - (A^A/C - (A/i n ) 2 ) - 2 (At< + #V n ) (S"^ - d M a n ) 



+^ (^„„ - d u A nfJi ) - d v A») 



(14) 
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Here A is the difference operator, denned in (1), which operates on the indices n, whiie 
At is defined as 

N-l 

A f 0n = ^ a k <l>n-k , (15) 

which satisfies the following equation: 

N-l N-l 



n=0 n=0 

The Lagrangian is invariant under transformations with the local parameters and 

hnnu > h n ^ u + d^ nv + di/^n^ , 

+ A f Cn • (17) 

Since the spectrum of A is given by (4), the Lagrangian describes a massless graviton 
(with mass m° = 0) and N — 1 massive gravitons with masses \m M \ (M = 1,2, • • • , N — 
1) in (4). We should point out that the massive gravitons always appear in pairs, 
sharing a common mass. As we have shown, the complex mass parameter m M (M = 
1, 2, • • • , N — 1) can be arbitrary chosen, just by properly choosing the coefficients in 
(1) to (6). As discussed in [5], the Lagrangian (14) can be regarded as corresponding 
to a Kaluza-Klein theory where the extra dimension is a lattice. 

As an example, we consider the two-brane RS model [6] (for a recent review, see [8]). 
In that model, the masses m M of the Kaluza-Klein modes are given by the solutions 
of the following equation: 



= (y' 2 (0) - |y 2 (0)) (j' 2 {z c ) + 



:J2 [Zc 



2(l + z c 

^)4i.(0))(^W + ^w) • (18) 

Here I is the length parameter of the five-dimensional AdS space and z c is given by the 
geodesic distance irr c between the two branes as 

z c = l(e^-l) . (19) 

The functions j 2 (z) and y 2 [z) are given in terms of Bessel functions as 

j 2 (z) = v^TT J 2 (m M z + l) , y 2 {z) = \fz~+l Y 2 (m M z + /) . (20) 

Specifically, we can consider models given by a first few solutions for m M of Eq. (18), 
in the different cases (8), (9), and (10). On the other hand, when m M is large, the 
Bessel functions behave as 

^zJ 2 {m M z) ~ \[^M cos (™ Mz ~ X*) > 

V-zY 2 (m M z) ~ y^Bin( TO ^-^) . (21) 
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Then Eq. (18) reduces to 



tan (m M l - K^j ~ tan (m M (z c + I) - ^tt ] . ( 22 ) 

which has the solution 

m M = — , (23) 

K being a positive integer. Motivated by (23), we may consider a N = 2N' + 1 graviton 
model with graviton masses given by 

w M_/^> M = 0,l,...iV', 

m M = AT' + l,7V' + 2,-.-iV-l = 2iV'. (24j 
Eq. (6) gives 

a = , 



2vr 



(2iV' + l)z c 



■ 2nN'n \ ■ 2irn 

1 _ e 2AT' + 1 e 2JV' + 1 



27rn 



„• 2 ( nnN' 

(-l) n 27T sm {-FT 



Nz c sin 



X — e *2jv'+i 



, (n / 0) . (25) 



Then iV plays the role of a cutoff of the Kaluza-Klein modes. 

In previous models [3, 5] of deconstruction, mainly nearest neighbour couplings 
between the sites of the lattice had been considered. Then, on imposing a periodic 
boundary condition, the lattice looks as a circle. Departing from this case, in the 
model considered here, we have introduced non-nearest-neighbour couplings among 
the sites. Then, a site links to a number of different ones in a rather complicated way. 
In this sense, the lattice in the present model is no more a simple circle but it looks 
like, say, a mesh or a net. Let us assume that the sites on the lattice correspond to 
points in a brane. If the codimension of the spacetime is one, the brane should be 
ordered, resembling the sheets of a book. One brane can only interact (directly) with 
the neighbouring two branes. However, if the spacetime is more complicated and/or the 
codimension is two or more, brane can (directly) interact with three or more branes, an 
interaction that will be perfectly described by our corresponding model. For example, 
a site on a tetrahedron connects with three neighbouring sites. In this way, the non- 
nearest couplings may adequately reflect the structure of the extra dimension. In this 
respect our model is very general and opens a number of interesting possibilities. 

3. The induced cosmological constant. We will here compute the one-loop ef- 
fective potential and as a consequence the induced cosmological constant for the ./V 
graviton model discussed in the previous section (for a recent review on Casimir energy 
calculations from extra dimensions, see [9]). We shall use zeta- function regularization 
[10], but instead, a different regularization scheme could work as well. First of all we 
compute the effective potential for a free scalar field with mass M since this corre- 
sponds to the contribution of each degree of freedom to the one-loop effective potential 
of our theory. 
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In the zeta-function regularization scheme, the one-loop contribution to the effective 
potential is given by 



V eff = -^C'(O) , (26) 

V being the volume of the manifold and £(s) the zeta function corresponding to the 
operator L = — V 2 + M 2 . In a flat manifold it can be easily computed by means of 

c(s) = J_ r dtt s-i Tve ~tL/^ = y£ r dtt s-d/2-i - t M* 

V ' r(s) Jo (4vr)<V 2 r( S ) J 

V M d T(s-d/2) ( M\~ 2s 



(4vr) d / 2 T(s) V » ) 



- 



(27) 



which is valid in arbitrary dimension d. The parameter has to be introduced for 
dimensional reasons and it will be fixed by renormalization. We should notice that 
the last expression provides the analytic continuation of the zeta function to the whole 
complex s— plane. One can see that in even dimensions there are a finite number of 
simple poles at the integer points s = 1, 2, d/2, while in odd dimensions there is an 
infinite number of simple poles at all the half-integer points s < d/2. In particular, in 
4-dimensions the zeta function has only two simple poles, and reads 

VM A f M\~ 2s 



32vr 2 (s - l)(s - 2) \ /i / 
By taking the derivative of Eq. (27) with respect to s, we obtain 



f(1) _ (-l) rf / 2 M c 



eff 2(d/2)!(47r) d / 2 
which is valid in even dimensions, and 



M 2 

In — k aid) 



(29) 



(1) _ MT(-d/2) 

V eff~ 2(47r) d/2 ' W 

valid in odd dimensions. Here a(d) = 1 + 1/2 + 1/3 + .. . + 2/d is a rational number 
larger than 1. 

It has to be noted that, in the odd dimensional case, the effective potential does 
not depend on /i and, owing to the presence of the Euler's gamma function, it is 
alternatively negative and positive. In particular, it is positive in five dimensions. 
The situation drastically changes in even dimensions. In fact in such case there is an 
explicit dependence on the scale parameter, and the effective potential can be positive 
or negative, this depending on the value of Mj [i. 

In the multi-graviton model which we are considering in this paper, the effective 
potential can be obtained by adding up several contributions of the kind (29) for d = 4. 
The masses are given by Eq. (24). In our model we have a massless graviton (which 
does not give contributions to the effective potential) and an even number of massive 
gravitons ( with the same mass). Then, in order to perform the computation, it is 
sufficient to consider only one-half of the whole mass spectrum, that is 

7T V N — 1 

z c z 
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In this way, we get 

N-1 . 

* M* f M 2 
V eff = VM + lO'E^fc In -| 



p=i 



iV-1 

2 _2,„4 / ^2^2 



p=l 

TW , 2- 10 iV(Af 2 -l)(3iV 2 -7)7r 2 / ^ 2 3 

7T 2 p 4 

64iJ 



JV-1 



p=i 



where Vr(/x) is the renormalized vacuum energy and the factor of 10 comes from the 
fact that there are two gravitons with the same mass, each of them having 5 scalar 
degrees of freedom. 

The effective potential is a physical observable and for this reason it has not to 
depend on the choice of the arbitrary parameter scale \i. This means that it has to 
satisfy to some renormalization condition (see e.g. [11]) 

„^-0, (33, 

from which it follows 

VM = Vr(hr) + 2 "° N{N2 - ff 3 * 2 " 7)7:2 In £ , (34) 



fj,R being the renormalization point which can be fixed by the condition Vr(ur) = 0. 
In this way one finally gets 

2- 10 iV(iV 2 - l)(3iV 2 - 7)vr 2 / vr 2 3 
Veff = 33 r^|"2 



c 



N-1 

2 J±J± 



p=l 

It interesting to note that, whatever the value of u R is (in practice it is fixed by some 
standard renormalization condition), due to the particular form of the mass spectrum, 
the expression (35) turns out to be positive if the number of gravitons is sufficiently 
large. This is a quite natural way to obtain a positive cosmological constant. One can 
observe that for large values of N, the effective potential goes as const • iV 5 (In N)/z*. 

Now let us consider (for comparison) the multigraviton model studied in Ref. [5]. 
Also in that case there is a massless graviton, which can be omitted, and N— 1 gravitons 
with masses 

M p = 2msin(^), p = 1, 2, N - 1 . (36) 

m being a positive constant. From (28) and (29) with d = 4, it directly follows that 
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and the corresponding effective potential, renormalized according to the procedure 
explained above, is 



where the factor 5 reflects again the scalar degrees of freedom of massive gravitons. 
The sign of the effective potential for this second example (Eq. (38)) can be positive or 
negative, depending on the value m/pR, but, contrary to what happens in our model, 
it is always negative, for the renormalization mass pr is large as compared with m (see 
[5]). The last expression has a nice behavior with N as N — > oo, namely ~ aN, being 
a = -0.041180. 

We should now remark that these results correspond to our re-calculation of the 
model in [5]. In fact, we noticed that a zeta function with a strange behavior was 
obtained there: apparently the starting point was a standard representation for C( s )i 
valid in 4 dimensions for Sfts > 2, but after analytic continuation a representation for 
C(s) valid only for 3is < 2 was obtained (Eqs. (16) and (22) in [5]). This was probably 
because of an unclear treatment of the massless graviton, which gives a divergent 
contribution in the original expression. In particular, owing to the existence of this 
zero mode, it seems that the step from Eqs. (19) to (20) in [5] missed the fact that the 
contribution of p = is special: all integer values of £ (and not only those which are 
multiple of N) contribute. This further yields a correction to Eqs. (22) and (23) in [5], 
which is exactly of the same kind as the last term there but with N = 1. Obviously, 
this new contribution is infinite for q = 1 (see the last term of Eq. (23) in [5]). This 
is why our Eq. (38) differs from the corresponding one in [5]), what is rendered clear 
when comparing the large iV behavior of both. 

4. Discussion. Let us now estimate the distance between the two branes in our multi- 
graviton theory obtained from a discretized Randall-Sundrum model in order to induce 
the observable cosmological constant. We may write V e ff as V e ff = Vq (N, z c , p) z~ 4 ; 
as the Hubble constant H is H ~ 10~ 33 eV in the present universe, the effective 
cosmological constant, A e //, should be A e ff ~ (lCT 33 eV) . Since the Planck scale is 



10 19 GeV~ 10 28 eV, one gets V eff ~ ^ ~ 10 -33x2+28x2 ( eV )4 = 1Q _i ( eV )4_ Then if 



length is 10 33 cm, the value of this z c corresponds to 10 3 to 10 2 cm. In the 
(continuum) model [12], if the number of extra dimensions is n, the length of the 

30 1 7 

extra dimension will be 10 « cm. Then, if z c corresponds to the length of the extra 
dimension, the value obtained corresponds to the n = 2 case. We note that the value 
we have obtained for z c is within the present limits of the measurements. 

As a consequence, the new multi-graviton theory we propose admits quite an arbi- 
trary mass spectrum. To be more specific, a model of this kind may be induced by a 
discretized brane-world. The induced cosmological constant is then defined by quan- 
tum effects. With a reasonable choice for the distance between the branes, the induced 
cosmological constant lies actually within the range of the admissible observational val- 
ues. Thus, a multi-graviton theory with non-nearest-neigbour couplings in the theory 
space can provide a useful model for dark energy in the current accelerating universe. 





(38) 
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The whole approach may open a new window in the realization of the deconstruc- 
tion program for QED/QCD like theories. Indeed, by similar adjustments of the mass 
spectrum one is able to obtain quite noticeable changes in the four-dimensional scalar 
effective potential. And such a dependence of the effective potential on the choice of 
mass spectrum allows us directly to put bounds to reasonable latticized boundary con- 
ditions from phenomenologically accepted estimations of the four-dimensional scalar 
potential. 
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